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Abstract 

The spectra of previously constructed auxiliary matrices for the six- vertex model 
at roots of unity are investigated for spin-chains of even and odd length. The two 
cases show remarkable differences. In particular, it is shown that for even roots of 
unity and an odd number of sites the eigenvalues contain two linear independent 
solutions to Baxter's TQ-equation corresponding to the Bethe ansatz equations 
above and below the equator. In contrast, one finds for even spin-chains only one 
linear independent solution and complete strings. The other main result is the 
proof of a previous conjecture on the degeneracies of the six-vertex model at roots 
of unity. The proof rests on the derivation of a functional equation for the auxiliary 
matrices which is closely related to a functional equation for the eight-vertex model 
conjectured by Fabricius and McCoy. 
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1 Introduction 



Baxter's TQ-equation [Q is one of the corner-stones of integrable systems and has been 
discussed in a variety of contexts. While it originated from the Bethe ansatz computa- 
tions for the six- vertex model, it provides conceptually a more general framework to solve 
the transfer matrix eigenvalue problem of integrable vertex models and paved the way 
for the solution of the eight-vertex model [HEl- (See also for recent developments.) 

In the TQ-equation the T symbolizes the transfer matrix of the integrable model 
at hand and Q is called the auxiliary matrix in terms of which the transfer matrix can 
be expressed. To be concrete, consider the six-vertex model on a square-lattice with 
periodic boundary conditions. Then the TQ-equation is a second order linear difference 
equation of the form 

T(z)Q(z) = ^zq-^Qizq 2 ) + <P(z)Q(zq- 2 ) (1) 

with z € C being the spectral variable and q the coupling (or crossing) parameter of 
the model. The coefficients <j) are given by some known scalar function, which in our 
convention specified below will be <j)(z) = (zq 2 — 1) M . Here M is the number of lattice 
columns. This equation is often treated on different levels. Some authors interpret it 
merely on the level of complex-valued functions, i.e. for the eigenvalues of the respective 
matrices only, not addressing the usually harder problem of the explicit construction 
of the matrix Q. Solving this construction problem, however, allows one to address 
the problem of determining the existence and number of possible solutions to the TQ- 
equation. This is particularly important for those cases which one often finds to be 
explicitly (or implicitly) excluded from the discussion, for instance when the parameter 
q is a root of unity or when the length of the spin-chain associated with the vertex model 
is odd. 

For these "special" cases the relation between the spectrum of explicitly constructed 
Q-matrices in El El El and the solutions of the TQ-equation will be our primary 
interest in this article. For our discussion it will be important to distinguish between 
the TQ-equation as a functional relation and an operator equation. 

1.1 The TQ-equation in terms of functions. 

In |lf)j Krichever et al pointed out that Baxter's TQ-equation naturally appears as an 
auxiliary linear problem in the context of the discretized Liouville equation, see Section 
4 therein. For the six-vertex model the analogue of this equation arises from the fusion 
hierarchy I12j. 

{zq 2 )T^ {z) - T< n+1 ) {z)T^ n - 1 ^ {zq 2 ) = (f>{z)^{zq 2n ) . (2) 

Here the T( n )'s should be thought of as the eigenvalues of the six- vertex fusion matrix 
of spin n/2 with = T being the transfer matrix and = 4> the quantum deter- 
minant |13j . We will give the definition of the six- vertex fusion hierarchy below. The 
general exposition given in J01 (which also addresses the case of higher rank and elliptic 
functions) concludes that there will be in general two-linear independent solutions, say 
Q ± , to the functional equation £Q) which satisfy the Wronskian 

Q + (zq 2 )Q-(z) - Q + {z)Q-{zq 2 ) = <j>(z) . (3) 
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From this Wronskian two types of "complementary" Bethe ansatz equations can be 
deduced [TU] . Notice that after equation (4.27) in ^U] the restriction to elliptic poly- 
nomials (j) of even degree is made which corresponds to the case of spin-chains of even 
length. While the focus in ^U] is mainly on the elliptic case, one wonders about the 
trigonometric limit and the implications for the six-vertex model. 

Pronko and Stroganov investigated this question of two potential linear independent 
solutions of the six- vertex TQ-equation in Jl] . Starting with the XXX spin-chain their 
discussion is generalized to the XXZ spin-chain excluding the root of unity case. They 
identify for spin-chains with an odd number of sites M one solution, say Q + , as the 
familiar solution from the Bethe ansatz, while the second solution Q~ is associated with 
solving the Bethe ansatz on the "wrong side of the equator" meaning that it incorporates 
M — n Bethe roots. Here < n < M/2 is the number of down-spins in the corresponding 
eigenstate of the transfer matrix or spin-chain Hamiltonian. For spin-chains with an even 
number of sites, however, there appears to be only one solution in the case of periodic 
boundary conditions. Baxter addresses this phenomenon on the basis of the coordinate 
Bethe ansatz and numerical computations in ^H] stating that the missing Bethe roots 
for the second solution Q~ have gone off to infinity. (See also the third paragraph after 
equation (4.43) in [EJ-) When a nonzero horizontal electric field is applied, numerical 
computations show both linear independent solutions exist for M even and odd |15j . 

The difference between spin-chains of odd and even length has been further under- 
lined in [T2|. For the special root of unity case q^ = 1 and spin-chains with an odd 
number of sites the two independent solutions of the six-vertex TQ-equation have been 
explicitly constructed ^H] starting from a conjecture on the particularly simple form of 
the groundstate [T71ITH1 ITU] . 

Our discussion in this article will show that the solutions obtained at q s = 1 in ^H] 
are in fact the only linear independent solutions in the spin S z = ±1/2 sector which 
possess the expected number of (M — l)/2 and (M + l)/2 Bethe roots, respectively. 
This is due to the large degeneracies in the spectrum of the transfer matrix connected 
with the loop algebra symmetry of the six- vertex model at roots of unity [2U] . A similar 
reduction in the number of solutions takes place at higher roots of unity as well. In other 
words, the number and nature of solutions to the TQ-equation (JJ) does also crucially 
depend on the value of the parameter q and not only on the length of the spin-chain. 

The various cases outlined and the necessary distinctions one has to make clearly 
show the importance of the explicit construction of auxiliary matrices in order to obtain 
better control over the different scenarios. However, Baxter's construction of auxiliary 
matrices for the six and eight-vertex model as presented in his book [21] applies only to 
spin-chains of even length (see the comment after equation (9.8.16)) and also at roots 
of unity can only be extended to a particular subset of cases. In contrast, the auxiliary 
matrices constructed in [7] for roots of unity q N = 1 and the ones in jf)] for the case 
of "generic" q do not have such limitations and apply for even as well as odd length 
of the chain. They are also of a simpler algebraic form. In this work we will explicitly 
relate those auxiliary matrices' spectrum to the two linear independent solutions of the 
TQ-equation by extending the discussion of [S] from the even to the odd case. 
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1.2 The TQ-equation in terms of operators. 



Starting with the papers \2'2\ \2'6\ on the Liouville model and subsequent papers on the 
six-vertex model [241 H3 [7] a new approach to construct auxiliary matrices has been 
developed which relies on representation theory. In this method one first solves the 
Yang-Baxter equation to obtain a matrix which commutes with the transfer matrix and 
afterwards derives the TQ-equation by investigating the decomposition of certain tensor 
products of representations. The novel feature [Hid in this context is the appearance of 
additional free parameters, collectively called p, in the auxiliary matrix, i.e. Q = Q(z;p). 
These free parameters shift in the operator solution of the TQ-equation [HUH) 



Thus, in this construction method one has to consider a generalized version of Baxter's 
TQ-equation at the level of operators. The free parameters are, for instance, necessary 
to break spin-reversal symmetry or to lift the degeneracies of the transfer matrix at roots 
of unity [Jj. They also contain information on the analytic structure of the eigenvalues 
of the auxiliary matrices [HI IH] ■ 

For the case when q is not a root of unity and twisted boundary conditions the 
spectra of of the auxiliary matrices constructed in [5] have been computed using the 
algebraic Bethe ansatz |5J. It was found that the eigenvalues decompose into two parts 
which are related by spin-reversal 



For special choices of the parameters p the functions Q are the eigenvalues of the 
lattice analogues of the Q-operators constructed by Bazhanov et al for the Liouville 
model (HI For twisted boundary conditions these two parts of the eigenvalue 

© can indeed be identified (up to some normalization constants) with the two linear 
independent solutions Q^ 1 of the TQ-equation discussed in the previous section, see 
[5] and the conclusions of this article. However, in general this is not true, since the 
two linear independent solutions Q^ 1 might not always exist, as for example in the 
case of periodic boundary conditions and spin-chains of even length. Notice also that 
Baxter's Q-operator for even spin-chains corresponds to only one of these eigenvalue 
parts, say Q + . This should make clear that the Q-operators in |)J[7] not only differ in 
the construction procedure but are quite different objects from Baxter's Q discussed in 



1.3 Outline and results of this article 

The outline of the article is as follows: 

Section 2. We introduce the fusion hierarchy of the six-vertex model and fix our con- 
ventions. The fusion matrices are defined such that they are polynomials of max- 
imal degree M in the spectral variable z. Here M is the length of the associated 
XXZ spin-chain. 

Section 3. We recall the definition of a particular subset of the auxiliary matrices at 
roots of unity constructed in |2j and discuss the general form of their spectra defin- 
ing the subparts of the eigenvalue decomposition (JSJ) . Preliminary results have 



T(z)Q(z;p) = ^(zq- 2 )Q(zq 2 ;p') + <j)(z)Q(zq- 2 ;p") . 



(4) 




(5) 



EH- 
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already been obtained in |Hj for the case of spin-chains with M even. Additional 
results for M even and odd are contained in [5]. In this article we will complete 
our investigation of the spectra by extending the discussion of [Hj to spin-chains 
with an odd number of sites. 

Section 4. The eigenvalues of the auxiliary matrices are in general polynomials in the 
spectral variable. In this section we determine when they have maximal degree 
and when they vanish at the origin. Both facts are related to the absence or 
occurrence of infinite Bethe roots. 

Section 5. We review from [S] the discussion of the TQ-equation in terms of the eigen- 
values of the auxiliary matrices. In particular, we discuss the transformation 
under spin-reversal and show how the second linear independent solution to the 
TQ-equation arises. 

Section 6. One of the main results of this article is the derivation of a functional 
equation, for M even and odd, that proves a previously formulated conjecture in 
(HI, see equations (18) and (19) therein. This result will enable us to determine 
the level of degeneracy of the eigenvalues of the transfer matrix at roots of unity 
and relate the two parts Q ± of the eigenvalues via an inversion formula. Moreover, 
this result can be considered as a six- vertex analogue of the eight-vertex functional 
equation conjectured by Fabricius and McCoy in We will comment on this 

in the conclusions. 

Section 7. The results for the case of odd spin-chains will show that for even roots of 
unity and periodic boundary conditions we obtain both linear independent solu- 
tions Q ± of the TQ-equation from the auxiliary matrices in [7j. This is a con- 
structive existence proof for these solutions. For odd roots of unity we will see 
that the eigenstates of the transfer matrix associated with (in the case that 
both solutions exist) correspond to zero eigenvalues of the auxiliary matrices. We 
in particular make contact with Stroganov's solutions for N = 3 ^B] an d show 
that they are the only ones with the expected number of Bethe roots (see also the 
appendix). 

Section 8. We state our conclusions and relate our results to the case when q is not 
a root of unity considered in |5j and the recent developments in the eight-vertex 
model IHIS1- 

2 The six-vertex fusion hierarchy 

We introduce the six- vertex model from a representation theoretic point of view. Denote 
by Wz : U q {sl2) — > EndC n+1 the spin n/2 evaluation representation of the quantum 
loop algebra, i.e. 




[n — m + l] q \m — 1} 
[m + l] g \m + 1} , 5 




TT 



i n) (/i)H 




TT 



(6) 
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with m = 0, 1, n. Define the fusion matrix of degree n + 1 by setting 

T (n + D {zq -n- 1 } = Tr • • • ^(Z^- 1 ) (7) 

where L( n+1 ) is the intertwiner with respect to the tensor product ir^ (g> ir^ , 

(0| L^ n+1 \w) |0) = ^vr^^^ 1 / 2 ) - ^( n )(g- hl / 2 ), 

(0|L( n+1 )H|l} = wqiq-q-^iq^TT^ih), 

(l\L^ n+1 \w)\0) = (q-q-^ie^iq-^ 2 ), 

(l|L (n+1) (u>) |1) = wq^ n \q- hl/2 ) - {q hl/2 ) . (8) 

Here the scalar products are taken in the second factor of the tensor product, i.e. the 
spin 1/2 representation. The fusion matrices satisfy the functional equation 

T^{z)T^{zq~ 2 ) = {zq 2 - l) M T( n+ V (zq' 2 ) + (z - 1) M T^^zq 2 ) . (9) 

The fusion hierarchy contains two special elements from which all others can be succes- 
sively generated, namely the six- vertex transfer matrix* T and the quantum determinant 
[TH] which are obtained via the identification 

T {2 \zq- 2 ) = T{z) and T^{z) = {zq 2 - 1) M . (10) 

In this manner the above functional equation may also serve as a defining relation for 
the fusion matrices. An alternative form of the fusion hierarchy is the one given in the 
introduction, see lf2"j). Both versions are equivalent. From the transfer matrix T we 
obtain as logarithmic derivative the XXZ spin-chain Hamiltonian, 



IT I ~X\ d 1 T ( Z ) 

Hxxz = -{q-q ) z— In 



dz {zq 

M 



(11) 



z=l 



I f q + q^ 1 1 

= ~2^2 } a m a m+l + ^1+1 H 2 (^m^rn+l ~ 1) \ ■ ( 12 ) 

m=l ^ J 

The well-known symmetries of the model are expressed in terms of the following com- 
mutators 

[ T ( m ){ z ),T^{w)] = [T^{z),S z ] = [rW (*),«*] = [tW(z),6] = , (13) 

where the respective operators are defined as 

M M M 

^ = -^C X = o*®-®o* = ]Ja x m , 6 = a z ® • • • ® a z = J] . (14) 

m=l m=l m=l 

These symmetries hold for spin-chains of even as well as odd length. 

*Note that our definition of the six- vertex transfer matrix differs from the one in 

El El 03 by an overaii 

factor, 

T{z)^ q ^T{z)/(zq 2 -1) M . 
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3 Auxiliary matrices at roots of unity 



In a series of papers |S1 1§] auxiliary matrices for the six- vertex model at roots of unity 
have been constructed. The construction procedure and the difference with Baxter's 
method have been discussed in [7] and we refer the reader to this work for details. 
In order to keep this paper self-contained we briefly give the definition of a special 
subvariety of the set of auxiliary matrices constructed in [7j . 

3.1 Definition 

Suppose q is a primitive root of unity of order N and set N' = N if the order is odd 
and N' = N/2 if it is even. Define the following ^'-dimensional nilpotent evaluation 
representation tt^ of the quantum loop algebra Ug^slz) [25| 126]. 



^{q hl )\n) = ^q- 2 "' 1 \n) , tt"(/i) \n) = \n + 1) , \N' - l) = , 



and 

K(q h0 )=7r»(q- h i), <(/„) = «rV( ei ), <(e ) = «; 7r"(/i), »,^C X . (16) 
Let the matrix 

L M = ( a ^ ? M J = a M ® cr+o-" + /? ® ct+ + 7 M cr~ + <^ <g> cr~o-+ . (17) 

be the intertwiner of the tensor product 7r^ (g)7r2=i of evaluation representations, explic- 
itly 

a M H = wq^^l 2 )-^-^' 2 ), f3^w)=wq(q-q- 1 )^(q h ^ 2 )^(f 1 ), 

7 M = (g-^^vr^exK^ 172 ), <VH = wqir»(q- hl/2 ) ~ ^(q hl/2 ) . (18) 

Define the auxiliary matrix in terms of these matrices as the trace of the following 
operator product, 

Q„(z) = TrL^ / (z/^ M _ 1 (^)---^ 1 (z/^) . (19) 

This matrix commutes by construction with the fusion matrices, 

[Q fi (w),T^(z)}=0, (20) 
and preserves two of the symmetries (fT3|) |H] , 

[Q fl (z),S z ] = [Q fl (z),e}=0. (21) 
Spin-reversal symmetry on the other hand is broken [7J |H| , 
XQ^Z, q)VK = Q^-i (z^ 2 , q) = Q^-i (zq 2 ^ 2 , q' 1 ) 1 = {-zq/ n) M Q^q' 2 . (22) 
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These relations hold for all M and allow one to determine the conjugate transpose of 
the auxiliary matrix [H], 

Q„(z, q)* = Qn(z, q- 1 ) 1 = Q-^zq- 2 , q) . (23) 

In addition, one derives from the following non-split exact sequence of evaluation rep- 
resentations TTw |7j 

-» ^ < ® 4 X) -> 7r ^»" 1 -» °> w = w 'i~ l = W "l = Z IV ( 24 ) 

the TQ-equation 

T(z)Q^z) = (z- l) M Q m {zq 2 ) + (zq 2 - l) M Q w -i(zq- 2 ) . (25) 

The proof can be found in |7j, here we will only review parts of the calculation of 
the spectrum of the auxiliary matrices given in [Hj and extend the results therein to 
spin-chains of odd length. 

3.2 The general form of the spectrum 

The starting point is the same as in jSJ: provided that the commutation relation 

[Q lx (z),Q u (w)} = 0, fx,v,z,weC (26) 

holds, the eigenvectors of Q^(z) are independent of the parameter \x as well as the spec- 
tral variable z. In order to prove (|26|) one has to explicitly construct the corresponding 
intertwiners of the tensor products iTw <S> k\ for all N' G N. As pointed out in [H] the 
necessary conditions for these intertwiners to exist are satisfied for all roots of unity. An 
explicit construction has been carried out for iV = 3 [HI and N = 6. Numerical checks 
have been performed for N' = 4, 5, 7. We shall take this as sufficient evidence for (|2H|) 
to hold true. 

There are two important implications of (j2(ifl . The first is that the auxiliary matrices 
are normal and hence diagonalizable, see (|23|) . The second consequence is that the 
eigenvalues of must be polynomial in the spectral variable z. Their most general 
form is therefore given byt 

Q»(z) = M^z n -P B {z)P il (z)P s (z N '^ 2N ') 

ra+ n- n$ 

= z n °° H(l - z/zi) ;Q(1 - *MM) 11(1 - z N '/a^)) . (27) 

i=l i=l i=l 

Note that we slightly differ in the notation from |Sj and have redefined the normalization 
constant by setting -Pb(O) = P^(0) = Ps(fy = 1- Besides these minor differences our 
definition of the various polynomials entering the eigenvalues is the same as in 8 . 

• The monomial factor in front of the eigenvalue is related to the occurrence of 
vanishing and infinite Bethe roots when the root of unity limit is taken in the 
deformation parameter q. 

^Throughout this paper we will denote eigenvalues and operators by the same symbol. 
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• The second factor Pb contains only roots Zi which do not depend on the free 
parameter [i and which will be identified with the finite Bethe roots at roots of 
unity. Moreover, we exclude from this set complete or exact strings, i.e. for every 
Zi there is at least one integer < I < N' such that Ziq 2e is not a root of Pb- 

• The third factor contains roots which do depend on the parameter \x. For even 
chains, M € 2N, this factor was identified in [Hj with the rescaled polynomial Pb, 

P^z) = P B (z^ 2 ) . 

Here we will find that this relation ceases to be valid for even roots of unity when 
M € 2N + 1. Again we exclude the possibility that the roots of P^ occur in strings. 

• Finally, the last factor P$ contains all roots which occur in strings. The zeroes 
dj may or may not depend on the parameter \x. Because all roots are sitting in a 
string the polynomial depends on z N rather than z. Note that we allow for the 
possibility n$ = 0. 

For later purposes let us decompose the eigenvalues of the auxiliary matrices into two 
parts similar as it has been done in [231 |Hj for q "generic" . Namely, we set 

n + 

Q + (z) = P B (z) = H(l-z/z i ) (28) 



and secondly, 



Q-(z) = lim Mfj, z n ~P^z)P s (z N ,fi 2N ). (29) 



This decomposition may seem arbitrary at the moment but it will become clear in our 
line of argument. Notice that we have eliminated the dependence of the auxiliary matrix 
on the free parameter [i in (|29|) . Moreover, it is worth stressing that Q are not always 
both identical with the two linear independent solutions Q 1 * 1 of the TQ-equation (the 
latter might even not exist), 



T(Z) ^ T^\zq- 2 ) = (zq* - 1) A V 52 ^^ + (* - 1) A V 5Z ^ 



(30) 



which satisfy the Wronskian 

q SZ Q + {zq 2 )Q~{z)-q- SZ Q+{z)Q~(zq 2 ) = {q SZ -q- SZ )(l- zq 2 ) M . (31) 

The additional phase factors q ±sz in comparison with the equation © of ^0] discussed in 
the introduction are due to different conventions. We shall set Q + (0) = Q~~ (0) = 1 and 
deg Q + = M — deg Q~ = M/2 — S z . While Q + is by definition the solution Q + above 
the equator, Q~ ^ Q~ in general. In particular we have included the normalization 
constant lim, n _ N > J\f u in the definition of Q~ , which in some cases can be zero as we 
will discuss below. 



8 



4 The degree of the eigenvalues and "infinite" Bethe roots 

We start our discussion with the first factor in the eigenvalue (|27|). the monomial z n °° 
which is related to the fact that some Bethe roots in the root of unity limit vanish 
or tend to infinity. Obviously, ^ if and only if the eigenvalue of the auxiliary 
matrices vanishes at the origin. Another obvious observation is that by construction of 
the auxiliary matrices it follows that 

deg Q tl = n OQ + n + +n-+ n s N' < M . (32) 

As it turns out this upper bound is assumed if and only if we have a vanishing monomial 
contribution, i.e. = 0. This can be deduced from the following relation for the 
auxiliary matrices |Hj 

Q^z) = (-zq/^Q^z^q-y (33) 
which implies that the coefficients in the power series expansion 

M 

QM = Qt^ m (34) 

m=0 

are related via 

qM = ( _ /u) -M g -M +2m (QW-™)y . (35 ) 

In particular, setting m = in the above identity we see that the eigenvalue of is of 
degree M whenever Q M (0) = Q$ ^ 0. At the same time this clearly prevents ^ 0. 
The coefficient Qjy can be easily calculated by noting that the building blocks L^(z) of 
the auxiliary matrix are lower triangular matrices at z = 0, 

LM(0) = -^{q- h ^ 2 ) ® a + a- - ^{q h ^ 2 ) ® a~a + + (q - q' 1 )^ {e x q- h ^ 2 ) ® <T . (36) 

Thus, the matrix Q^O) in quantum space is diagonal with its diagonal elements given 
by 

N'-l 

q„(o) = n M Tr q - h ^ = n M (^r £ ^ ■ ( 37 ) 

For the interpretation of this result we distinguish the following cases: 

1. When q N ' = 1, i.e. for primitive roots of unity of odd order, the degree of the 
polynomial will only be equal to M in the commensurate sectors 2S Z = OmodiV. 
At the same time this means that there are no infinite Bethe roots, that is = 0. 
Consequently, we have 

Q M (0) =M» = (—) M N q s y sz , q N ' = 1, 2S Z = OmodiV . (38) 

2. When q N ' = -1 we have to distinguish between M even and odd. Let M be 
even then S z takes integer values only and we obtain deg = M if and only if 
2S Z = mod N. Again we find in these spin-sectors the normalization constant 

Q M (0) =Af (t = N'q sz fj, sz , q N ' = — 1, M G 2N, 2S Z = OmodiV . (39) 
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3. For q N = — 1 and M odd, however, the total spin eigenvalue will be constrained 
to the set 2S Z £ 2Z + 1, preventing the existence of a monomial factor. That is, in 
this case we always have deg = M and = 0. The normalization constant is 
therefore 

N' S z —N'S Z 

Q,(0)=M, = -q N ' S ^ SZ q ~ q s , , q N '=~l, M£2N+1. (40) 

5 The TQ-equation 

The most important property of the auxiliary matrix is the solution of the following 
functional equation with the six-vertex transfer matrix which has been proved to hold 
for M even and odd 0, 

T(z)Q tl (z) = (z- l) M Q m {zq 2 ) + (zq 2 - 1) M Q w -i(zq~ 2 ) . (41) 

From this functional equation and the fact that does not depend on the free pa- 
rameter fi we infer similar to the case M even considered in (Hj that 

Wi(/j,) = Wiji 2 and cii(fJ>) = ai({j, 2N ) 

implying the following form for the eigenvalues of the transfer matrix, 

m = ^ ( , _ 91^1 + ^ _ W) . (42) 

Here the ratios of the normalization constants can only depend on q from which we 
deduce 

^ = A . (43 ) 

The zeroes of the polynomial Q + are fixed through the "Bethe ansatz" equations, 

= (1 - z iq 2 ) M q- s Q + (z iq - 2 ) + (1 - z t ) M q s Q + (z t q 2 ), i = 1, ...,n+, (44) 

with 

q s ;= *m q 2n x = "o^f*- . (45) 

We will argue below that this phase factor is determined by the total spin of the eigen- 
state and the number of Bethe roots which tend to zero in the root of unity limit. 



5.1 Spin reversal and Bethe roots "beyond the equator" 

We are now exploring the role of the polynomial factor P^. Following the same line 
of argument as in (Bj we act with the spin reversal operator from both sides on the 
TQ-equation employing the transformation law [7| |S] 

mQ II (z)K = Q, l -i(z f i- 2 ) 
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of the auxiliary matrix. Replacing afterwards 1 we obtain the equation 

T{z)Q^ 2 ) = (z- l) M Q m -,{zn 2 ) + (zq 2 - l) M Q m {zn 2 ) . (46) 

Since we employed the spin-reversal operator we refer to this identity as the TQ-equation 
"beyond the equator" . The corresponding expression in terms of eigenvalues is deduced 
to be 

[ ) ~ [Z lj P,{z^) W M [zq L) P M (V) 

= (* - + (V - 1) V • (47) 

Q (z) Q (z) 

Notice that in comparison with 1)42(1 not only the phase factors q s have been inverted but 
that the polynomial Q~ will in general have a different degree than Q + , i.e. n_ ^ n + . 
From 1)46)1 we now obtain the "Bethe ansatz equations beyond the equator" , 

0=(l-w i q 2 ) M q s Q-(w i q- 2 ) + (l-w i ) M q- s Q-(w i q 2 ), i = l,...,n_. (48) 

This second solution to the TQ-equation, which is related by spin-reversal to Q + , will 
not always be linear independent as we shall see below. The parameter s entering the 
phase factors in the eigenvalue expressions of the transfer matrices changes sign, since we 
already saw for particular cases, see (|38j). (j3HJ), (0OJ), that it is related it to the total spin 
of the corresponding eigenvectors through the normalization constant. For the general 
case its value can be determined by making contact with the fusion hierarchy. 

Using the functional equation @ for the fusion matrices presented in the introduc- 
tion, the above results for the transfer matrix are extended to all fusion matrices via 
induction. A straightforward calculation yields 



r <„> (2) _ 5± <.«>. Q±W W) vj . (49) 



In the paper [S] it has been argued, using the algebraic Bethe ansatz when q is not a 
root of unity, that the parameter s can be identified with 

s = 2n + S z modN' . (50) 

Here no denotes the number of Bethe roots which vanish in the root of unity limit 
q N - 1. 



6 A functional equation relating Q + and Q 

The final step in the analysis of the spectrum of the auxiliary matrices rests on the 
following functional equation, which has been proved for iV = 3 in [H], 

Q^ 2 q 2 )Q v {z) = (^ 2 -l) M g^(Vg 2 )+g jV ' M Q Mi/? -^ +1 (V9 2 )^ JV '- 1 )(^ 2 ) • (51) 

This equation is a direct consequence of the following decomposition of the tensor prod- 
uct ttw (g) of evaluation representations, 

- < (2) < - 7T { f- 2) - (52) 
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where 

2 1 i II —N'+l II -N'+l i N'+l 

w = \xvq , \i = (ivq, w = fj,q, \i = [ivq , w = \iq , z = vq T . 

(53) 

For the moment assume the functional equation 1)51(1 to hold, the derivation of (|51jl and 
(f52*|) is given in the appendix. Let us insert the explicit form of the eigenvalue into (|51|) . 
We find 

Q^z^q^Qyjz) ^ 

^ z ^Pn(z)p . rwh ft (1 " ^V 2Af/ M(M 2jV '))(i - ^M(^)) 

M„ uq ^iWW Ul (1 _ z N> (J 2N>/ a .( JJI 2N> l ,2N>) 

M ^f' +1 gN'MrpiN'-l)^ + (^2 _ 1} M 

Here we have exploited the fact that the zeroes of the factors P^, P u and P^n'+i, 
P^vq only depend on [i 2 ,v 2 and n 2 v 2 q 2 , respectively. The possible complete string 
contribution P$ contains the various parameters only to the power 2N' . Notice that 
the last line of the above equation is independent of the free parameters fj,, v (the ratio 
of the normalization constants only depends on q as pointed out earlier). This implies 
that the zeroes ai(n 2N ) are either independent of \i altogether or only incorporate it 
as a multiplicative factor, i.e. one has (exactly as in the case N = 3 proved in |HJ) the 
alternative 

a i ( f j, 2N ')=a i or a^') = a { f N ' ■ (54) 

Hence, there are at most 2 ns possible eigenvalues of the auxiliary matrix in a degenerate 
eigenspace of the transfer matrix with fixed tZqo and n±. This proves part of the second 
conjecture made in [S], see equation (18) with aj = (zf) N . From the outcome on the 
zeroes of the complete string contribution one deduces that the factor originating from 
the Ps polynomials simplifies, 

(1 - z N 'ii 2N '/ ai (ii 2N 'v 2N ') ~ S{ ^~ L >- 

Hence, we can rewrite the functional equation in terms of Q as follows, 

~- 2n -Q + {z)Q-{zq 2 ) 



Q^(zfi 2 q 2 )Q u (z) Njt 

Q/iuq(zfi 2 q 2 ) Af -N'+l 



q K'M T ( N, -V(zq 2 ) + (zq 2 - l) M 



^g-N' + l 



Here we have set v = q N in the normalization constants without loss of generality. 
Invoking now the earlier stated form of the fusion matrices (|49|) this identity becomes 

q- 2n -Q + (z)Q-(zq 2 ) = -*fB— q -s Q +( z )Q-( zq 2 ) = 



Nftq-N'+l Af M -N' + l 



+ a N M±N s Q ± (z)Q ± (za 2 ) V - { — - J + (zq 2 - 1" 



M 
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Setting fj, — > q N and solving the last expression for we obtain the result 

N ' n -2is( vr M -\\M 

fW-WME jg^, . (50) 

Hence, the two solutions of the TQ-equation are related to each other by an inversion 
formula (provided that Q~ ^ 0, see the discussion below). This identity is the six- vertex 
analogue of the functional equation conjectured by Fabricius and McCoy for Baxter's 
1972 auxiliary matrix of the eight- vertex model [HIE]- (We will comment further on this 
in the conclusion.) Here we have proved this functional equation in the six- vertex limit 
for all roots of unity and spin-chains of even as well as odd length. Let us investigate 
the difference between the solutions depending on the cases when M is even or odd. 



Decomposition of the eigenvalue for M even. For spin-chains of even length 
M = 2m and at q N = ±1 the sum in (|29|). which at first sight appears to be a rational 
function, simplifies due to the Bethe ansatz equations (|44(1 to a polynomial, i.e. we have 
for any contour C encircling the point ziq~ 2 ^ , 



V ' -Ust„„2t 1 \M 



J- I y g n - *) dz 

2 ™ / h Q+(zq 2i )Q+(zq 2i - 2 ) 



+ -^n „ X , t— r = 0, < £ < N 



Q + (z iq - 2 ) n(i-^Ai) Q + (ziq 2 )U( 1 -*i/*u 

As a consequence we obtain in general the identifications [Sj 

lim P^{z) = Q + {z) (57) 



-JV' 



and 



We therefore conclude that Q~ ^ Q~ in this case. In fact, one finds numerically that the 
solution Qr with M/2 + S z Bethe roots does not exist. The other important conclusion 
is that the result (|58|) together with Q54|) enables us to read off the degeneracy of the 
transfer matrix eigenvalue corresponding to Q + . According to (|54[) each zero cij in the 
complete string contribution P$ is either independent of fi or is multiplied by a factor 
p? N showing that there are 2 ns possible eigenvalues of the auxiliary matrices each 
corresponding to a vector in the degenerate eigenspace of the transfer matrix. This is 
in agreement with the observation [20] that only spin-1/2 representations occur in the 
tensor products describing the finite-dimensional representations of the loop algebra. 
Obviously, if ns = the eigenvalue of the transfer matrix is non-degenerate up to 
spin-reversal symmetry. 
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Decomposition of the eigenvalue for M odd. For M = 2m + 1 odd the above 
simplification of the eigenvalue in general also holds true when N odd with the possible 
exception that the part (|58|) of the eigenvalue completely vanishes, we will discuss this 
case below. For N even, however, the polynomials P^ and Pb = Q + differ. Since we 
have now s € and q N = — 1 the rational function (|58|) has poles at z = z; t with 
non- vanishing residue, 

1 f N> a- 2ls (-\ yn 2t \ M 

2ni J ^ Q+izq^Q+izq^- 2 ) 
C(zi) £ " 1 

q- 2s (z t q 2 - 1) M q^Hzj - l) M 



Q + (^? 2 )ii(i-^Ai) Q + (^- 2 )n(i-^i) 

Therefore, the factor Q + (z) in front of the sum is needed to cancel these poles and the 
above factorization (|57|). (f55|) does not take place. As a consequence complete strings 
are absent and there are no additional degeneracies other than spin-reversal symmetry. 
As argued earlier in Section 4 infinite Bethe roots are absent as well and we must have, 

N ' n -2tS z ( vn 2l _ i\A/ 



^Q+(zqM)Q+(zq2 



with the degrees of the polynomials Q obeying 



degQ" = M-degQ+ = y + S z . (60) 

Thus, we obtain a very different picture depending on the length of the spin-chain being 
odd or even. 



The quantum Wronskian. The difference between the two situations of even and 
odd spin-chains is highlighted further by introducing the analogue of for the two 
different parts of the auxiliary matrix eigenvalues. This corresponds to the "quantum 
Wronskian" in |23| . First note that using Q56|) we easily obtain 

q ns Q + (zq 2n )Q-(z) - q- ns Q + (z)Q-(zq 2n ) = (q N ' s - q- N ' s )T^(z) . (61) 

Upon specializing to n = 1 this relation simplifies to 

q s Q + (zq 2 )Q-(z) - q- s Q + (z)Q-(zq 2 ) = (q N ' s - q- N ' s ){zq 2 - l) M . (62) 

Notice that the right hand of the above equation always vanishes except for odd spin- 
chains and even roots of unity. This signals the linear dependence between for M 
even and q N = 1, M odd as described above, compare with (|57|1 and (|58J1 . 

For M odd and q N = — 1 the quantum Wronskian is non-zero and we can identify 

N'S Z _ -N'S Z 

Q+ = Q+ and Q~ = M^.* Q~ = - g , _ g _ s , QT • (63) 
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Thus, via an explicit construction of diagonalizable Q-operators we have shown existence 
of the solutions above and below the equator. Notice that the Wronskian implies the 
Bethe ansatz equations (|44|). Namely, we have for each zero Z\ of Q + that 

^ Q+(V) _ q N>S*_ q -N'S* _ 

5 1\M ' V' 34 '' 



* {zrf-l)**- Q-( Zi ) * { Zi -1) M 

An analogous relation holds for the zeroes W{ of Q~ leading to the Bethe ansatz equations 
beyond the equator (]48[). Note that (|31|). respectively (|62[). contains more information 
than each copy of the Bethe ansatz equations by itself, as it relates the zeroes Z{ and W{ 
through the following sum rules for each < m < M, 

A/A n S z -2£ p + p - _ n -S z ~2k p + p - 

Here we have in light of Q63JI identified the zeroes of Q* 1 with {zj} and {wi} and intro- 
duced the elementary symmetric polynomials 

e+ = e fe (^ 1 ,...,^) and = e k (w^ 1 , wj}_ n+ ) . (66) 

Numerically it is by far more feasible to solve this set of M equations, which is quadratic 
in the M variables {e^~ } U {e^}, rather than the original ra+ = M/2 — S z Bethe ansatz 
equations (|44j) which are of order M in the n+ variables {zi}. We verified for iV' = 3, 5 
up to spin-chains of length M = 11 that the number of solutions of the equations ()65|) 
matches the dimension of the respective eigenspaces of the transfer matrix. 



7 Zero eigenvalues at TV and M odd 

An additional aspect in which the cases of even and odd spin-chains differ is the oc- 
currence of zero eigenvalues of the auxiliary matrices. That the auxiliary matrices can 
have indeed a non-trivial kernel for odd M has already been remarked upon in [7j where 
it was noted for the simple case of the M = 3 spin-chain. As we will see it is closely 
connected with the inversion formula (|56[) which follows from the functional equation 
(foT|) and the two independent solutions Q ± of the TQ-equation. 



Eigenstates with a maximum number of Bethe roots. Suppose q is an odd root 
of unity then, as discussed above, the rational function (jf)<S[) becomes a polynomial in 
z and Pfi{z) = Q + (zfi~ 2 ). The only exception to this scenario is the case when the 
corresponding eigenstate of the transfer matrix is a singlet. According to our previous 
discussion, we therefore must have 2 ns = 1, i.e. the corresponding eigenvalue of the 
auxiliary matrix cannot contain complete strings. In the absence of infinite Bethe roots 
we now argue that this implies (|58|) vanishes. This can be deduced in several ways. 
Suppose the two linear independent solutions Q^ 1 to the TQ-equation exist. Then they 
have to obey the quantum Wronskian ()31|) . Solving the latter for Q~ we obtain 

_ 2S > Q+(z)Q~{zqZ) _ 2y (l-^ 

^ [Z) ~ q Q+(zq2) +[ 9 ' Q+(V) ' 1 j 
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Iteration of this formula yields after A-steps 

ar(z) = g- ms ° Q g } + (* 2SZ - m + (z) ^ 9 ( q } 



Q+(zq 2N ) j^Q + (zq 2l )Q + (zq 2t ~ 2 ) 
which upon invoking the root of unity condition q N = q N = 1 gives 

This fact together with the identification Q + = Q + and (|56[) implies the vanishing of the 
corresponding eigenvalue of the transfer matrix. This does not mean that the second 
linear independent solution Q~ does not exist, it simply states that the normalization 
constant A/^ = i in the definition (|29|) of Q~ ^ Q~ is zero. Note also that (|68|) applies 
to non-degenerate states only, which decrease in number as M S> N due to the loop 
algebra symmetry at roots of unity |20| . 



M 


3 


5 


7 


9 


N = 3 


1/3 


1/10 


1/35 


1/126 


N = 5 


3/3 


8/10 


21/35 


55/126 


N = 7 


3/3 


10/10 


33/35 


108/126 



Table 1. Shown are the number of "maximal" solutions to the Bethe equations (i.e. 
n_|_ = M/2 — S z Bethe roots above and n_ = M — n+ below the equator with S z = 1/2) over 
the dimension of the spin-1/2 sector. The deformation parameter is chosen to be q = exp(27ri / N) . 
For 7V=3,5 it has been checked that the number of "maximal Bethe states" matches the dimen- 
sion of the kernel of the auxiliary matrix. 

In fact, for N = 3 in the spin-sector S z = ±1/2 there is only one state with the 
expected number of Bethe roots above and below the equator, the groundstate. This is 
in agreement with the results in [T^]. However, our starting point is different from the 
one in ^H]- Instead of making a conjecture on the explicit form of the groundstate of the 
six-vertex model, we simply start from the assumption that there exists an eigenstate 
with m = (M — 1)/2 Bethe roots in the spin S z = 1/2 sector. According to (|32j) complete 
strings cannot be present and thus must be a constant. But because of (f6*8"|) with 
Q + = Q + this constant is vanishing and we have the difference equation, 

(1 - z) M Q + (zq 2 ) + q-\l - zq 2 ) M Q + (zq' 2 ) + q- 2 (l - zq- 2 ) M Q + {z) = . (69) 

As our conventions differ from Stroganov's we review his calculation in the appendix and 
show that (|69|) has a unique solution which can be expressed in terms of hypergeometric 
functions. The same holds true for the second linear independent solution Q~ which 
has m + 1 roots. 

In the case of general G 2N+1 similar difference equations follow. Let M = 2m+l, 
then the eigenvalues of singlet states with = ns = in the spin S z = 1/2 sector 
with N > 3 satisfy 

N N-2 

£ = o, f N (z) = (i - z) M n Q ± (^ 21 ) 

t=\ 1=1 
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implying the following sum rules in terms of the elementary symmetric polynomials (|66l) 
in the m Bethe roots above and the m + 1 Bethe roots below the equator, 

o= E ( M k ) E nY 2i H- (to) 

k+l=n ^ ' fciH hfcjv_ 2 =' J ' — 1 

Here the integer n takes all values in the range 

< n = — ^ 1 mod N < N M ^ 1 ± 1 (71) 
and the different summation variables run over the intervals, 

< k < M, < kj < Tl . (72) 

In general the set of equations ()70|) is of order N — 2 and only for AT = 3 becomes linear 
in the variables {e^ }, where the equations are particularly simple to solve. Nevertheless, 
these sum rules are still an advantage over the Bethe ansatz equations which are of order 
M. 



8 Conclusions 

Let us summarize the new results obtained for the six-vertex model at roots of unity. 
First of all the discussion has been extended from even to odd spin-chains exploiting 
that the construction procedure for the auxiliary matrices in [7] does not have the 
same limitations as the one in Baxter's book |21j . This allowed us to reveal the major 
differences in the spectrum of the six-vertex model at roots of unity between these two 
cases. 

1. When the length of the spin-chain is even there are degeneracies in the spectrum 
of the transfer matrix for all roots of unity. These degeneracies are reflected in 
the spectrum of the auxiliary matrices through the occurrence of "complete string 
factors", see (|58|) in the text. The number of these strings, i.e. the degree n$ of the 
polynomial Ps in (|2*T|). determines the degeneracy of the corresponding eigenspace 
of the transfer matrix to be 2™ s . This is in accordance with the observations 
made in [20] . In order to arrive at this result we made use of the crucial functional 
equation (|51|) which severely restricts the dependence of the complete string factors 
on the free parameter /i entering the definition of the auxiliary matrix (|19j) . In 
addition, we employed (|5T|) to prove the identity ((S5|) which states that the string 
factors are determined (up to their dependence on the aforementioned parameter 
/j.) by the solution to the Bethe ansatz equations and the number of infinite Bethe 
roots which fix the eigenvalue of the transfer matrix. These results had previously 
been proved for N = 3 only and conjectured to hold true for iV > 3 [S]. Moreover, 
we deduced that the Bethe roots appear twice in the eigenvalue of the auxiliary 
matrices, once in the factor Pg = Q + and once multiplied by the factor fi 2 in 
the factor of the eigenvalue (|27j) . A second linear independent solution to the 
TQ-equation was not found. 
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2. For spin-chains with an odd number of sites the novel feature was the appearance 
of such a second linear independent solution to the TQ-equation below the equator. 

For primitive roots of unity of odd order this second solution does not exist 
for all eigenstates of the transfer matrix, but only for those which are singlets 
and have no infinite roots. For these eigenstates of the transfer matrix we have 
shown that due to the functional equation (|51|) the corresponding eigenvalues of 
the auxiliary matrices must vanish, see equation (|68[) in the text. The number of 
these states, i.e. the dimension of the kernel of the auxiliary matrix, will become 
smaller as the length of the spin-chain starts to exceed the order of the root of 
unity, i.e. M S> N. (For N = 3 we in particular saw that there is only one 
such singlet state for all odd M and it corresponds to Stroganov's solutions of 
the TQ-equation for the groundstate of the XXZ spin-chain at A = —1/2 |16j.) 
This decrease in number can be understood in terms of the loop algebra symmetry 
|20| of the six-vertex transfer matrix. As the length of the spin-chain M grows 
more and more of the transfer matrix' eigenstates organize into larger and larger 
multiplets spanning the irreducible representations. Similar to the case of even 
spin-chains these degenerate states within the multiplets give rise to complete 
strings in the eigenvalues of the auxiliary matrices with the same formula 2™ s 
yielding the multiplicity of the transfer matrix eigenvalue. 

For primitive roots of unity of even order the solution below the equator 
always exists, here, however, the eigenvalues of the transfer matrix do not vanish. 
We showed the absence of infinite Bethe roots as well as complete strings, leaving 
at most a double degeneracy in the spectrum of the transfer matrix due to spin- 
reversal symmetry. The latter is broken by the auxiliary matrices and we used 
this fact to identify Pb = Q + and P^ = Q~ with the solutions to the TQ-equation 
above and below the equator, respectively. The explicit construction of the Q- 
matrices in |7J, see the definition (|19|) in this article, guarantees therefore the 
existence of these two solutions, a fact implicitly assumed in Jl] for the case of 
"generic q" . What is still lacking at the moment is an understanding of the physical 
significance behind the existence of two linear independent solutions opposed to 
the case when there is only one. We hope to address this question in a future 
publication. 

In this article we have focussed on the case when q is a root of unity to discuss the 
spectra of the auxiliary matrices constructed in [7J. But as mentioned in the introduction 
analogous Q-operators have also been constructed when q is not a root of unity Their 
spectra together with the resolution of certain convergence problems originating from an 
infinite-dimensional auxiliary space have been discussed in [2] . As explained therein one 
in general needs to impose quasi-periodic boundary conditions on the lattice in order to 
obtain a well-defined auxiliary matrix. For instance, in the critical regime \q\ =1 the 
twist parameter A has to be of modulus smaller than one, |A| < 1, to guarantee absolute 
convergence (U] • This twist parameter enters the definition of the fusion matrices ((7J) as 

T^ n+l \zq- n - 1 ) = Tr\* (n) ^® l L^ l \zq n+l ) ■ ■ ■ (zq^ 1 ) (73) 
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and modifies the Wronskian ()31j) in the following manner, 

A-V z Q+(V)Q-(2) - \q~ SZ Q + (z)Q-(zq 2 ) = (A" 1 / 2 - A(T 5 *)(1 - zq 2 ) M . (74) 

Numerical computations show that both solutions Q^ 1 exist for M even and odd. Pro- 
ceeding similar as we did for the case of odd roots of unity we can iteratively solve this 
equation for, say Q + , to obtain 



Q~(z) = lim A 



2n 

n^oo 



2n\ 



Q+(zq 2n ) 



,_ () Q+(^+2 )Q+(zg 2£ ) • 

In the limit n — > oo the first term on the right hand side tends to zero as |A| < 1 and 
|g| = 1. The above expression then matches the results obtained for the eigenvalues of 
the Q-operator from the algebraic Bethe ansatz, see equations (75-78) in 

Q<(z;r , n ) = {-l) M rf\l-\ 2 q - 2SZ )Q + (z ri )Q-(z) . 

Here, up to some trivial normalization factors, we can identify oc Q^ 1 . Thus, similar 
as in the root-of-unity case the spectrum of the auxiliary matrices (HI E] when q is not 
a root of unity contains both solutions to the TQ-equation, the one above and the one 
below the equator. Moreover, these two solutions are related by the analogue of the 
formula (|56|) where the summation extends now over an infinite interval. 

In the text we commented on a similarity between the relation (j56|) and a eight- 
vertex functional equation conjectured by Fabricius and McCoy |3J |Hj for Baxter's 1972 
auxiliary matrix Q$ v 2 at coupling values rj = mK/N' (see equation (3.10) in [I] or 
(3.1) in 0), 

e -inuM/2K Qgv ( u _ iK ,^ = 

in f h M (u-(2t + l)K/N>) 

W ] to Q ^ u ~ 2 ^/N')Q 8v (u - (21 + l)K/N>) ' 

Here h(u) = 9^(0)8i(ttu/2K)6^(ttu/2K) in terms of Jacobi's theta-functions with mod- 
ular parameter p = exp(—nK'/K). We have made the replacements L — > N' (the order 
of the root of unity) and N — » M (the length of the spin-chain) in the notation of [IJ [5] . 
For general N' the explicit form of the constant A is as yet unknown. This functional 
equation has been proved for the free fermion case when M is even and numerically ver- 
ified for coupling values corresponding to roots of unity of order three, see the comment 
after (3.1) in 0. 

If one identifies Qs v {u — iK') with Q~{e u q~ 1 ) and Qg v (u) with Q + (e u q~ v ) in the 
six-vertex limit the similarity becomes apparent. This is further supported by the ob- 
servation that the transformation u — > u — iK' corresponds to spin-reversal in the 
eight-vertex Boltzmann weights^. It is tempting to speculate on further identities such 

*The author is thankful to Barry McCoy for discussions on this point. 
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as the elliptic analogue of the relation 1)62(1 for instance. While these identifications can 
be conjectured and numerically investigated on the level of eigenvalues, the analogous 
construction of the eight-vertex Q-matrices corresponding to (|19|) . which would allow 
one to prove existence, is a more complicated problem. 
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A Derivation of the functional equation 

As mentioned earlier the proof of the functional equation (|51|) employs representation 
theory and is deduced from the following non-split exact sequence describing the de- 
composition of the tensor product iTw <g> tt" , 

- 4' ^ < ® < ^ 4' ® *f ~ 2) - . (76) 

The various parameters appearing in the representations are not all independent but 
satisfy the relations 

2 / / // —N'+l a N'+l i N'+l 

w = fiuq , fi = fJLvq, w = fiq, \i = fiuq , w = fiq , z = vq T . 

(77) 

Notice that we have set the second evaluation parameter in the tensor product Ww ® ^\ 
equal to one. The general case -Kw (g) -n v u is obtained by simply replacing w — > urn, w' — > 
w'u and w" — > w"u, z' — > z'u. The line of argument is analogous to the one applied in 
[Z] to derive the TQ-equation via (|24(1 and the proof of (|51j) in [H] for N = 3, whence we 
will be rather brief in presenting the various steps of the proof. 

The inclusion. First we determine the subrepresentation 7r^, contained in the tensor 
product TTw®^i wn en w is tuned to the special value given in (|52l) . This will yield the first 
term on the right hand side of (|51j) . The corresponding inclusion map i : ir^, =— > Tiw <8> 7r^ 
is determined by identifying the lowest weight vectors in both representations, 

7#3|0>^|0)®|0>e<®<. (78) 

The remaining relations for the rest of the vectors in the included subrepresentation is 
obtained by successive action of the quantum group generators via the formula 

|0> ^ «®<)A(x)|0)®|0) . 

The above formula suffices to compute the various parameters. For instance, the pa- 
rameter fi' in (|77[) is computed from the action of the Cartan element x = q hl . The 
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evaluation parameters w, w' are deduced as follows. First act with fi on the lowest 
weight vector to obtain 

<' (/i) |0) = \l)^> « ® <)A(A) |0) ® |0) = j/g |1> ® 10) + |0) ® |1> . 
Alternatively, one obtains via the generator eo, 

(eo) 10) = w' «, (8 <)A(e ) [0) |0> = w |1> ® |0) + M |0) ® |1> . 

Comparing both results leads to the stated values of w, w' in (|77|). The scalar coefficient 
in front of the first term on the right hand side of the functional equation (|51jl is obtained 
from the identity 

L^ 3 (zfiq 2 )L^(z/u)(t ® 1) = (V - 1)0 ® l)ir«(znq/v) (79) 
which can be easily verified by acting on the lowest weight vector. 

The projection. In order to compute the second term of the functional equation (|51j) 
we need to determine the representation in the quotient space ir^ ® tt^/vt^, which is 
projected out by the map r. One proceeds in an analogous manner. The projection 
map t is fixed by identifying this time the highest weight vector in both representation 
spaces, 

< ® < 9 |iV - 1) ® |AT' - 1) ^ |iV' - 1) ® | AT' - 2) G tt£, tt^'' 2) . (80) 

Again the value of //" can be inferred from the action of the Cartan element x = q hl via 
the formula 

« ® <)A(x) |iV' - 1) ® | AT' - 1) ^ «" ® 7r^'- 2) )A(^) \N' - l) ® |AT' - 2> . 

Setting x = fifo and x = eoei one obtains the desired evaluation parameters u/' and z' 
detailed in Q77|) . For instance, from the left hand side of the above equation one obtains 

(q - T^K, ® <)A(/i/ ) |iV - 1) ® |iV' - 1) = 



m + m 1 - m 2 



M 1?2 + i/ + z.- 1 - - v~\ 2 \ \N' - 1) ® | AT' - 1) 



while the right hand side is computed to 

(q - ^ 1 ) 2 «" ® ^'" 2) )A(/i/o) |iV' - 1) ® |AT' - 2> = 
^q-N'+i + ^-l q N'-l _ ^q-N'-i _ 



| AT' - 1) ® | AT' - 2) 



+ ( g -r'm--2] g|JV< _ 1)8K _ 2) 

Matching the coefficients in both results yields the stated expressions for the parameters. 
In the case of the quotient projection there is only a trivial additional scalar factor as 
we have the equality 

(r ® (W)J&(*A0 = q N 'L^- N ' + \z M - N ' +1 ML^- 2 \z q N ' +1 )(r ® 1) . (81) 

Again this is most easily calculated by acting with both sides of the equation on the 
highest weight vector. This completes the proof of the functional equation. 
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B Stroganov's solution revisited 



As explained in the text the assumption that there exists an eigenstate with m = 
(M — l)/2 Bethe roots in the spin S z = 1/2 sector implies via (|52|). (|55|) and (jHEJ with 
Q + = Q + that we have the difference equation, 

(1 - z) M Q + (zq 2 ) + q-\l - zq 2 ) M Q + (zq- 2 ) + q~ 2 (l - zq- 2 ) M Q + {z) = . (82) 

Expanding 

3m+l 

(l-z) M Q + {zq 2 ) = l+ J2 c n z n 

n=l 

we infer that the difference equation implies 

c n = if n = 2 mod 3 . 

The remaining coefficients can be determined from the fact that (1 — z) M Q + (zq 2 ) has 
an M-fold zero at z = 1. Applying the method of Lagrange interpolating polynomials, 
similar as it has been done in |27l I16j . one finds the ratios 

C3n+3 _ (n - m)(n - m - 1/3) c 3w+4 _ (n - m)(n + 1/3 - m) 

c 3n ~ (n + l)(n + 2/3) ' ,U< c 3ri +i ~ (n + l)(n + 4/3) 
together with 

c = 1 and ci = -(4/3) m /(2/3) m . 

Here (x) m is the Pochhammer symbol. From the ratios of the coefficients we infer that 
there is a unique solution which can be expressed in terms of hypergeometric functions 

(1 - z) M Q + (zq 2 ) = 2 F 1 (-m, -m - ±, §; z 3 ) - ^ z 2 iq(I - m, -m, |; z 3 ) . (83) 

(3)774 

Note that this solution is not simply obtained by multiplying Stroganov's solution (11) 
in ^H] with an exponential factor. This is due to the fact that we solved the difference 
equation (|69|) in terms of polynomials which are regular at origin, while Stroganov's 
solution applies to Laurent series. 

Second solution. Besides the solution for the Bethe polynomial we just obtained, 
there is a second solution "beyond the equator" as it possesses m + 1 roots. Set 

3m+2 

(i-,) m q-(V) = i+ E c ^ n 

n=l 

then it obeys the difference equation with S z = —1/2, i.e. 

(1 - z) M Q-(zq 2 ) + q(l - zq 2 ) M Q-(zq- 2 ) + q 2 (l - zq- 2 ) M Q~ (z) = . 

This implies for the coefficients 

c' n = if n = 1 mod 3 . 

As before the solution to this set of equations can be expressed in terms of hypergeo- 
metric functions, 

(l-z) M Q~(zq 2 )= 2 F 1 (-m,-m-l 1 r ,z 3 )-^z 2 2 F i q-m,-m,l;z 3 ) . (84) 

( 3 )m 
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Groundstate eigenvalue. From the difference equation it follows that the eigenvalue 
of the transfer matrix is given by 

T(z)P B = q±kz - 1) M ^^ + <^(V " l) ^^ = (zq-* ~ 1) M 

which matches the conjecture [171 1181 ITU] employed in [T^]. The corresponding ground- 
state eigenvalue of the Hamiltonian is given by Hxxz = —M . 
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